Quantum many-body dynamics in optomechanical arrays 
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We study the nonlinear driven dissipative quantum dynamics of an array of optomechanical sys- 
tems. At each site of such an array, a localized mechanical mode interacts with a laser-driven cavity 
mode via radiation pressure, and both photons and phonons can hop between neighboring sites. 
The competition between coherent interaction and dissipation gives rise to a rich phase diagram 
characterizing the optical and mechanical many-body states. For weak intercellular coupling, the 
mechanical motion at different sites is incoherent due to the influence of quantum noise. When 
increasing the coupling strength, however, we observe a phase transition towards a regime of phase- 
coherent mechanical oscillations. This transition and the phase diagram of the system are studied 
using a Gutzwiller ansatz for the dynamics of the driven-dissipative system. 



Introduction. - Recent experimental progress has 
brought optomechanical systems into the quantum 
regime: A single mechanical mode interacting with a 
laser-driven cavity field has been cooled to the ground 
state [TJ [2]. Several of these setups, in particular op- 
tomechanical crystals, offer the potential to be scaled up 
to form optomechanical arrays. Applications of such ar- 
rays for quantum information processing [3[ |4] have been 
proposed lately. Given these developments, one is led 
to explore quantum many-body effects in optomechani- 
cal arrays. In this work, we analyze the nonlinear photon 
and phonon dynamics in a homogeneous two-dimensional 
optomechanical array. In contrast to earlier works [3H6], 
here we study the array's quantum dynamics beyond a 
quadratic Hamiltonian. To tackle the non-equilibrium 
many-body problem of this nonlinear dissipative system, 
we employ a mean-field approach for the collective dy- 
namics. First, we discuss photon statistics in the array, 
in particular how the photon blockade effect [7] is al- 
tered in the presence of intercellular coupling. The main 
part of the article focusses on the transition of the collec- 
tive mechanical motion from an incoherent state (due to 
quantum noise) to an ordered state with phase-coherent 
mechanical oscillations. For these dynamics, the dissipa- 
tive effects induced by the optical modes play a crucial 
role. On the one hand, they allow the mechanical modes 
to settle into self-induced oscillations [5UT5] once the op- 
tomechanical amplification rate exceeds the intrinsic me- 
chanical damping, see Fig. [ljb) . On the other hand, 
the fundamental quantum noise (e.g. cavity shot noise) 
diffuses the mechanical phases and prevents the mechan- 
ical modes from synchronizing. This interplay leads to 
an elaborate phase diagram characterizing the transition. 
To gain further insight, we develop a semiclassical model 
describing the coupling of the mechanical phases and the 
influence of quantum noise. 

While true long-range order is prohibited for a 
two-dimensional system with continuous symmetry, 
a Beresinskii-Kosterlitz-Thouless transition towards a 
state with quasi-long range order is possible. The ordered 
mechanical phase thus resembles the superfluid phase in 




Figure 1. (a) Optomechanical array with localized mechani- 
cal (bj) and laser-driven optical modes (dj) at each site. The 
optical and mechanical coupling between neighboring sites is 
set by J and K, respectively, (b) Onset of self-induced oscil- 
lations for an isolated mechanical mode as a function of laser 
driving strength (schematic). The classical dynamics (black 
solid line) show a bifurcation. Quantum fluctuations blur the 
transition (dashed blue line) and generate a mechanical state 
whose phase is completely undetermined, see also Fig. [3jb) . 

two dimensional cold atomic gases [16] or Josephson junc- 
tion arrays [17]. Notably, optomechanical arrays com- 
bine the tunability of optical systems with the robustness 
and durability of an integrated solid-state device. Other 
driven dissipative systems that have been studied with 
regard to phase transitions recently include cold atomic 
gases [18-23 , nonlinear cavity arrays [24] [25] and opti- 
cal fibres [26 . In a very recent work and along the lines 
of [18] , the preparation of long-range order for photonic 
modes was proposed using the linear dissipative effects 
in an optomechanical array [6 . Our work adds the novel 
aspect of a mechanical phase transition to the studies of 
driven dissipative many-body systems. 

Model - We study the collective quantum dynamics of 
a two-dimensional homogeneous array of optomechanical 
cells (Fig. [I]). Each of these cells consists of a mechanical 
mode and a laser driven optical mode that interact via 
the radiation pressure coupling at a rate go (h = 1): 

Homj = -Aa]a J - + ^6]6 J --^o(^+^)4^+ a ^(4 + ^')- 
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The mechanical mode (bj) is characterized by a frequency 
Q. The cavity mode (dj) is transformed into the frame 
rotating at the laser frequency (A = cji aser — cj cav ) and 
driven at the rate a^. In the most general case, both 
photons and phonons can tunnel between neighboring 
sites (ij) at rates J/z and K/z, where z denotes the 
coordination number. The full Hamiltonian of the array 
is given by H = J2j H om ,j + #int, with 

To bring this many-body problem into a treatable form, 
we apply the Gutzwiller ansatz A\Aj « (A\)Aj + 
A\{Aj) - (A\)(Aj) to Eq. {2). The accuracy of this ap- 
proximation improves if the number of neighboring sites z 
increases. For identical cells, the index j can be dropped 
and the Hamiltonian reduces to a sum of independent 
contributions, each of which is described by 

Hmf = H om - J (a! (a) + a(a+)) - K(&(b) + 6<6+».(3) 

Hence, a Lindblad master equation for the single cell den- 
sity matrix p, dp/dt = —i[H m t, p] + nV[a)p + TV[b]p can 
be employed. The Lindblad terms V[A)p = ApA^ — 
\A^ Ap — T^pA^A take into account photon decay at a 
rate n and mechanical dissipation (here assumed due to 
a zero temperature bath) at a rate V. 

Photon statistics. - Recently, it was shown that the ef- 
fect of photon blockade [7 can appear in a single optome- 
chanical cell: The interaction with the mechanical mode 
induces an effective nonlinearity for the photon field of 
strength g^/Q. [3 [27]. Hence, the presence of a single 
photon can hinder other photons from entering the cav- 
ity. To observe this effect, the nonlinearity must be com- 
parable to the cavity decay rate, i.e. g^/^t > and the 
laser drive weak (az, <C k) [7} [28]. 

To study nonclassical effects in the photon statistics, 
we analyze the steady-state photon correlation function 
g^(r) = (at(t)at(t + r)a(t + r)a(t))/(a(t)ta(t)) 2 [29] at 
equal times (r = 0), with g^ 2 \0) = 1 for a coherent state, 
and g^ 2 \0) < 1 (> 1) indicating anti-bunching (bunch- 
ing). Here (Fig. [2|, we probe the influence of the collec- 
tive dynamics by varying the optical coupling strength J, 
while keeping the mechanical coupling K zero for clarity. 
We note that, when increasing J, the optical resonance 
effectively shifts: A — >> A + J . To keep the photon 
number fixed while increasing J, the detuning has to be 
adapted [30]. In this setting, we observe that the inter- 
action between the cells suppresses anti-bunching (inset 
of Fig. [5]). Photon blockade is lost if the intercellular 
coupling becomes larger than the effective nonlinearity, 
2 J > g^/ft. Above this value, the photon statistics shows 
bunching, and ultimately reaches Poissonian statistics for 
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Figure 2. Loss of photon blockade for increasing optical cou- 
pling in an array of optomechanical cavities. The equal time 
photon correlation function shows anti-bunching (g^ 2 \0) < 1) 
and bunching (g^ 2 \0) > 1) as a function of detuning A and 
optical coupling strength J. The smallest values of g^ 2 \0) 
are found for a detuning A = -go /SI. When increasing the 
coupling J while keeping the intracavity photon number con- 
stant, i.e. along the dashed line, photon blockade is lost (inset, 
g^ 2 \0) as black solid line). For a smaller driving power (inset, 
blue solid line, oil — 5 • 10 -5 ft), anti-bunching is more pro- 
nounced, and the behavior is comparable to that of a nonlin- 
ear cavity (inset, dashed line). The hatched area in the main 
figure outlines a region where a transition towards coherent 
mechanical oscillations has set in (see main text and further 
figures). ^ = 0.3 ft, a L = 0.65k, g = 0.5ft, r = 0.074ft. 

large couplings. Similar physics has recently been ana- 
lyzed for coupled qubit-cavity arrays, [30]. For very large 
coupling strengths, though, the density plot of Fig. [2] re- 
veals signs of the collective mechanical motion (hatched 
area). There we observe the correlation function to os- 
cillate (at the mechanical frequency) and to show strong 
bunching. We will now investigate this effect. 

Collective mechanical quantum effects. - To describe 
the collective mechanical motion of the array, we focus 
on the case of purely mechanical intercellular coupling 
(K > 0, J = 0) for simplicity. Note, though, that the 
effect is also observable for optically coupled arrays, as 
discussed above. 

As our main result, Figs. [3^ a) and |4ja) show the 
sharp transition between incoherent self-oscillations and 
a phase-coherent collective mechanical state as a function 
of both laser detuning A and coupling strength K: In the 
regime of self-induced oscillations, the phonon number 
(b^b) reaches a finite value. Yet, the expectation value 
(b) remains small and constant in time. When increas- 
ing the intercellular coupling, though, (b) suddenly starts 
oscillating: 

(b)(t)=b + re- in ^ t . (4) 
Our more detailed analysis (see below) indicates that 
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Figure 3. Transition from the incoherent to the synchronized 
(coherent) phase: (a) Mechanical coherence r, measuring the 
amplitude of the order parameter in (&)(£) = b-\-re~ ineiit , as a 
function of laser detuning A and mechanical coupling K. At 
weak coupling, the self-oscillations are incoherent, r = 0, due 
to quantum noise. When increasing the coupling strength, the 
systems shows a sharp transition towards the ordered regime, 
where the mechanical oscillations are phase- coherent, r > 0. 
A cut for fixed detuning (dashed line) is shown in Fig. [4^a). 
(b,c) Modulus of the density matrix elements (in Fock space) 
and Wigner density of the collective mechanical state in the 
incoherent (b) and the coherent regime (c) (as marked in (a)). 
go = K = 0.3 O, oll 1.1 «, T 0.074 CI 

this transition results from the competition between the 
fundamental quantum noise of the system and the ten- 
dency of phase locking between the coupled nonlinear 
oscillators. Below threshold, the quantum noise from the 
phonon bath and the optical fields diffuses the mechan- 
ical phases at different sites and drives the mechanical 
motion into an incoherent mixed state. The reduced 
density matrix p^ is predominantly occupied on the 
diagonal, see Fig. [3^d), and the Wigner distribution, 
W(x,p) = ^ f^oo( x ~ y\p^\x J ry)e 2wy / h dy, has a ring- 
like shape, reflecting the fact that the phase of the motion 
is undetermined p2J|3T], see Fig. [3^b). Above threshold, 
the mechanical motion at different sites becomes phase 
locked, and the coherence parameter r (Eq. Q) reaches a 
finite value. Within our mean- field ansatz, r grows with 
a critical exponent 1/2 close to the threshold, see Fig. 
(4^a). The emergence of coherence also becomes a_ppar- 
ent from the off-diagonal elements of p^ (Fig. p[e)). 
The corresponding Wigner function assumes the shape 
of a coherent state with a definite phase oscillating in 
phase space, see Fig. j3^c) . Thus, this transition spon- 
taneously breaks the time-translation symmetry. In a 
two-dimensional implementation, true long range order 
is excluded, but the coherence between different sites is 
expected to decay as a power law with distance. We 
also note that this transition is the quantum mechanical 
analogon of classical synchronization, which was studied 
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Figure 4. (a) The mechanical coherence r as a function of cou- 
pling strength K for fixed detuning A = 0.511. The dashed 
line shows the optical readout of coherence, i.e. the oscillat- 
ing component of the photon number (at a), which is directly 
accessible in experiment. The loss of synchronization for even 
larger coupling strengths can be attributed to enhanced phase 
fluctuations and a nonmonotonic behavior of the effective 
phase coupling rate between the phases, as discussed in the 
main text, (b) The phonon number (b^b) shows maxima at 
the resonance and at the sideband (A — Ao ~ Cl). (c) The 
diffusion constant for the mechanical phase, D^, for an un- 
coupled (K = 0, solid line) and coupled array (K = 0.1 Q, 
dash-dotted line). Other parameters as in Fig. [3] 



for optomechanical systems in [32j [33]. An important 
difference is, though, that the classical nonlinear dynam- 
ics were analyzed for an inhomogeneous (with disordered 
mechanical frequencies) system in the absence of noise 
[32j [33] , while in our case disorder is only introduced via 
fundamental quantum noise. 

The laser detuning determines both the strength of 
the self-oscillations and the influence of the cavity shot 
noise on the mechanical motion. As we will show be- 
low, even the coherent coupling between the mechanical 
phases (ultimately leading to synchronization) is tunable 
via the laser frequency. As a result, the synchronization 
threshold depends non-trivially on the detuning param- 
eter A, see Fig. [3|a). The onset of self-induced oscilla- 
tions happens continuously and already at a red-detuned 
laser frequency, A < 0, due to the effect of the static 
mechanical shift and the shift of the optical spectrum 
by Ao = — g^/Cl. Close to the onset, the mechanical 
phases are very susceptible to quantum noise preventing 
synchronization. The diffusion rate D^, defined as the 
linewidth of the correlator (b(t)&(0)) ~ e -^^+D v )t^ 
shows a maximum (Fig. [4^c), solid line). For finite cou- 
pling strengths, the diffusion is enhanced, most strikingly 
at the mechanical sideband. As a result, the synchro- 
nization threshold shows a minimum between the onset 
of self-oscillations and the sideband. From extended sim- 
ulations, we find that this behavior is generic for systems 
in the resolved sideband regime (Cl > k). 

Semi- classical description. - In order to gain further 
insight into the coupling and decoherence mechanisms, 
we analyze a semi-classical model for coupled Hopf oscil- 
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lators with phases cpj and amplitudes Aj (in units of the 
mechanical ground state width) . We follow [32 , but add 
noise terms: 



(fj — — \l — \a\ cos (fj 



K 



K 



i COS (fi COS ifj 



Aj = j(Aj -A) + ^2 A i cos W sin <Pj + (5) 

The radiation pressure forces the mechanical motion onto 
a limit cycle with steady state amplitude A and ampli- 
tude decay rate 7, and also modifies the mechanical os- 
cillation frequency via an amplitude dependent optical 
spring effect: Q(A) = Q — 2go(\a\ 2 cos (p)r/A, where (...)t 
denotes the average over one mechanical period. The 
fluctuating noise forces £^ and £4 comprise the effects of 
the phonon bath and cavity shot noise [34]. The Hopf 
equations ([5| can be derived directly from the Langevin 
equations of the optomechanical system, and a descrip- 
tion of this kind has proven good qualitative agreement 
with the full quantum dynamics for a single optomechan- 
ical cell [I2j [34] . The complexity of equations ([5| can be 
reduced by integrating out the amplitude dynamics [32] . 
In addition, we apply a mean- field approximation: We 
replace e lipj for neighboring cells by (e lipj ) = Re 1 ^ and 
e t2(pj ^Zffjj = R 2 e l ^2^ w here (...) denotes the average 
over all sites [35 . The resulting equation describes the 
coupling of the mechanical phase on a single site, (/?, to a 
mean field \I/: 



ip = —Q(A) + KRcos(^ -ip) + KxRsm^ - <p) 

+ K 2 R 2 sin(2^ - 2(p) + K 2 RR 2 sin(^ 2 -9-(p)+^ 

(6) 

where we defined the coupling rates K 2 = K 2 /2j and 
K\ = KdQ/j — K 2 and a fluctuating noise force 
which is characterized by a diffusion constant [34] , 
to be discussed below. The differential dVL = A^\ a =a 
accounts for the variation of the mechanical frequency 
due to fluctuations in the amplitude. It turns out to be 
decisive, since it determines the strength and the sign of 
the (most important) coupling rate K\. 

The phase equation ^ has the form of the mean-field 
Kuramoto equation [35] in the presence of noise, with 
additional coupling terms proportional to K and K 2 . In 
the incoherent regime, the order parameter R (R 2 ) is zero 
and the phase fluctuates freely. In the coherent regime, 
i.e. for < R, R 2 < 1, on the other hand, the different 
coupling terms force the phase tp towards a fixed phase 
relation with Note that within the mean- field approx- 
imation, only a configuration with ip = \I/ is stable. This 
state can be achieved via the coupling terms proportional 
to K\ and K 2l which provide a restoring force ~ K\^ 2 R 



acting on cp. The cosine term only renormalizes the os- 
cillation frequency. This statement can be clarified by a 
linear stability analysis following the approach of Stro- 
gatz et al. [36] for the generic Kuramoto model. At the 
threshold of synchronization, the time evolution of the 
small increment SR on top of the incoherent background 
(R = 0) is governed by SR = (K x - 2D ip )SR + 0(5R 2 ), 
where the coupling terms K and K 2 enter only as higher 
order contributions. As a result, the incoherent phase 
becomes unstable for 



K x = 2D, 



(7) 



defining the threshold of the transition. Moreover, if K\ 
becomes negative, no stable phase synchronization is pos- 
sible. This situation arises if dVt < 0, or for large intercel- 
lular coupling rates K > 2dQ, as observed in Fig. Qa). 

The total phase diffusion rate is given for vanishing 
coupling (K = 0) by = ±(D^ + 6 -^D A ) [34], 
where the diffusion rates D^a correspond to the noise 
terms |^,a(£) of Eq. ^ and are defined by 2D^ y A = 
J dr(^^A{t + t)<^,a(^))t- For the explicit expressions 
we refer to [34], and discuss the main qualitative feature 
here: Close to the threshold of self-oscillations, i.e. for 
small amplitudes and small amplitude damping 7, the 
phase is highly susceptible to noise. Moreover, while a 
large frequency differential SQ and a small value of 7 
increase the strength of the effective coupling rate Ki, 
this effect is overwhelmed by the quadratically enhanced 
diffusion rate. Thus, the semiclassical analysis affirms 
our observations of the quantum dynamics of the sys- 
tem. We note, though, that the semiclassical equations 
neglect higher-order correlations and can not fully cap- 
ture the quantum dynamics in the regime of strong cou- 
pling, go > k. 

Experimental prospects. - We note that observation of 
the mechanical phase transition does not require single 



photon strong coupling 



> 



k): The quantum fluctu- 



ations of the light field will dominate over thermal fluc- 
tuations as long as gQ\a\ 2 / k > ksT/Q. This is essen- 
tially the condition for ground-state cooling, which has 
been achieved utilizing high-Q mechanical resonators and 
cryogenic cooling [TJ [2] . In contrast, we note that the 
photon-blockade effect (Fig. [2| requires low tempera- 
tures T and g$ > Vtn, or at least, in a slightly modified 
setup [37, 38 , go > n. While still challenging, optome- 
chanical systems are approaching this regime [39] . 

The flexibility to design two dimensional arrays has 
been demonstrated for optomechanical crystals in [40] , 
and sufficiently strong optical and mechanical hopping 
rates are feasible, as shown for comparable structures 
in [32]. The strongest challenge would likely be the si- 
multaneous optical driving of many cells, although sim- 
ilar physics may be observed for many cells coupling to 
one extended optical mode [5 , 32 , 33 (thereby effectively 
realizing global coupling). We expect the transition to 



5 



be robust against disorder [32]. One could also study 
the formation of vortices and other topological defects 
induced by engineered irregularities and periodic vari- 
ations, and explore various different lattice structures. 
Thus, optomechanical arrays provide a novel, integrated 
and tunable platform for studies of quantum many body 
effects. 
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